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This paper presents an analysis of observational data on the p-mode spectrum 
of the star aCen B, and a comparison with theoretical computations of the 
stochastic excitation and damping of the modes. We find that at frequencies 
> 4500 yuHz, the model damping rates appear to be too weak to explain the 
observed shape of the power spectral density of aCen B. The conclusion rests 
on the assumption that most of the disagreement is due to problems modelling 
the damping rates, not the excitation rates, of the modes. This assumption is 
supported by a parallel analysis of BiSON Sun-as-a-star data, for which it is 
possible to use analysis of very long timeseries to place tight constraints on the 
assumption. The BiSON analysis shows that there is a similar high-frequency 
disagreement between theory and observation in the Sun. 

We demonstrate that by using suitable comparisons of theory and observation 
it is possible to make inference on the dependence of the p-mode linewidths on 
frequency, without directly measuring those linewidths, even though the a Ccn B 
dataset is only a few nights long. Use of independent measures from a previous 
study of the a Cen B linewidths in two parts of its spectrum also allows us to 
cahbrate our linewidth estimates for the star. The resulting cahbrated linewidth 
curve looks similar to a frequency-scaled version of its solar cousin, with the 
scaling factor equal to the ratio of the respective acoustic cut-off frequencies of 
the two stars. The ratio of the frequencies at which the onset of high-frequency 
problems is seen in both stars is also given approximately by the same scaling 
factor. 

Subject headings: stars: oscillations - stars: activity - Sun: activity - Sun: 
helioseismology - data analysis 

1. Introduction 

Stars like the Sun, which have sub-surface convection zones, display a rich spectrum of 
acoustic (p-mode) oscillations. The oscillations are stochastically excited and damped by 
the convection, and this gives rise to an extremely rich spectrum of modes. Measurement of 
the amplitudes and damping rates of the p modes therefore gives important information for 
constraining theories of convection in stellar interiors. 

Chaphn et al. (2007) used model computations of the excitation and damping of Sun- 
hke oscillations to make predictions of the p-mode spectra of a selection of stars on the lower 
main sequence. The computations revealed an extremely interesting feature in the predicted 
appearance of the spectra of those model stars that had effective temperatures cooler than 
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about 5400 K: the modelled power spectral density of the modes showed two maxima, at 
different frequencies. Chaplin et al. found a pronounced dip in mode power between the 
maxima when the computations were made for young stars; since the maxima are well 
separated in frequency, the predicted spectra took on a "double humped" appearance. In 
older main-sequence stars the dip was found to be much less pronounced, and instead the 
spectra showed a broad plateau of power. 

The Kl V main-sequence star a Cen B (HR5460) is a suitable candidate to test the 
predictions, since its effective temperature lies on the cool side of the 5400-K boundary given 
by the model computations, and data are available on its p-mode spectrum from observations 
made by Kjeldsen et al. (2005). As we shall demonstrate, the model computations predict 
a broad plateau of power for a Cen B; however, this broad plateau is not seen in the p- 
mode data of Kjeldsen et al. This paper reports on attempts to try and understand this 
disagreement between theory and observation. 

The observational data for our study are the aforementioned Doppler velocity observa- 
tions of a Cen B, made by Kjeldsen et al. The theoretical predictions we use are pulsation 
computations of the stochastic excitation rates and the damping rates of the radial p modes 
of the star. By making judicious comparisons of the theoretical computations and the ob- 
servational data we show it is actually possible to make inference on the p-mode linewidths 
of a Cen B, without directly measuring those linewidths, even though the Kjeldsen et al. 
observations span only a few nights. The conclusions drawn do rest on one important as- 
sumption: that most of the disagreement that is observed between theoretical predictions 
and observations of the radial-mode amplitudes of a Cen B resides in errors in the model 
damping rates, and not errors in the model excitation rates. We provide evidence in support 
of this rationalization from a similar comparison of theoretical and observational data of the 
p-mode spectrum of the "Sun as a star". The Sun-as-a-star data come from observations 
made by the ground-based Birmingham Solar-Oscillations Network (BiSON) (Chaplin et al. 
1996) in Doppler velocity. 

The layout of our paper is as follows. We begin in Section [2] with a description of 
the excitation and damping-rate computations. We then proceed in Section [3] to compare 
theoretical predictions and observations of the mode amplitudes, for both a Cen B and the 
Sun. Then, in Section HJ we show how discrepancies between the theoretical and observed 
amplitudes may be explained largely in terms of errors in the theoretically computed damping 
rates. We also demonstrate how inference may be made on the mode linewidths from use 
of the observed amplitudes and the theoretical excitation and damping computations. We 
finish in Section [5] with a brief discussion of the main points of the paper. 



-4- 



2. Predictions from analytical model computations 

The stellar equilibrium and pulsation computations that we performed are as described 
by Balmforth (1992a) and Houdek et. al (1999). These computations gave estimates of the 
powers and damping rates of the radial p modes of a Cen B. 

The model computations required four general input parameters to specify the stellar 
model: the mass, M; radius, R] effective temperature. Teg; and chemical composition. We 
took values for a Cen B of: M = 0.934 ± 0.007 M©; R = 0.863 ± 0.005 R©; and T^s = 
5288 ± 38 K (see Yildiz 2007, and references therein). The composition was fixed at X = 0.7 
and Z = 0.02, as for the model computations performed in Chaplin et al. (2005, 2007). 
While this composition differs from estimates given in the literature for the a Cen B system 
(of Z ^ 0.025; again, see Yildiz 2007), changes in composition at this level have only a 
second-order impact on the model calculated p-mode excitation and damping rates, and any 
changes that might be relevant here are smaller than the observational uncertainties (e.g., 
see Fig. 17 of Houdek et al. 1999). 

2.1. Pulsation computations 

Computation of the excitation rates and damping rates of the p modes demands a 
description of how the pulsations interact with the convection. This requires computation 
of the turbulent fluxes associated with the convective heat and momentum transport. These 
turbulent fluxes are obtained from a nonlocal, time-dependent generalization of the mixing- 
length formulation of Cough (1977a, b), with a mixing length calibrated to the Sun. In this 
generalization there is a parameter, $, which specifles the shape of the convective eddies. 
Then there are two parameters, a and b, which control respectively the spatial coherence 
of the ensemble of eddies contributing to the turbulent fluxes of heat and momentum and 
the degree to which the turbulent fluxes are coupled to the local stratiflcation. These two 
parameters control the degree of 'non- locality' of convection; low values imply highly nonlocal 
solutions, and in the limit a, 6 oo the system of equations formally reduces to the local 
formulation (except near the boundaries of the convection zone, where the local equations 
are singular). Cough (1977a) has suggested theoretical estimates for their values, but it is 
likely that the standard mixing-length assumption of assigning a unique scale to turbulent 
eddies at any given location causes too much smoothing; accordingly, somewhat larger values 
probably yield more realistic results. In this paper we therefore adopt two sets of values for 
the non-local parameters which provide reasonable results for the Sun, and other Sun-like 
stars: a? = h"^ = 600; and a? = h"^ = 300. We comment further on the impact of this choice 
in Section [231 below. 
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2.1.1. Envelope and pulsation models 

Both the envelope and pulsation computations assumed the three-dimensional Edding- 
ton approximation to radiative transfer (Unno & Spiegel 1966). The integration was carried 
out inwards, starting at an optical depth of r=10^^ and ending at a radius fraction r/R = 0.2. 
The opacities were obtained from the OPAL tables (Iglesias & Rogers 1996), supplemented 
at low temperature by tables from Kurucz (1991). The equation of state included a detailed 
treatment of the ionization of C, N, and O, and a treatment of the first ionization of the next 
seven most abundant elements (Christensen-Dalsgaard 1982), as well as 'pressure ionization' 
by the method of Eggleton, Faulkner & Flannery (1973); electrons were treated with rela- 
tivistic Fermi-Dirac statistics. Perfectly reflective mechanical and thermal outer boundary 
conditions in the pulsation computation were applied at the temperature minimum in the 
manner of Baker & Kippenhahn (1965). At the base of the model envelope the conditions of 
adiabaticity and vanishing displacement were imposed. Only radial p modes were considered. 

2.1.2. Stochastic excitation model 

The amplitudes of stochastically excited oscillations are obtained in the manner of Chap- 
lin et al. (2005). The procedure is based on the formulation by Balmforth (1992b) and Gol- 
dreich & Keeley (1977) but includes a consistent treatment of the anisotropy parameter $ 
of the turbulent velocity field (for details see the discussion in Chaplin et al. 2005). For the 
anisotropy parameter we adopt the value $ = 1.13, a value that was also considered by Chap- 
lin et al. The excitation model assumes a description in which the largest, energy-bearing 
eddies are described by the mixing-length approach. The small-scale convection is modelled 
by a turbulence spectrum, for which we adopt the Kolmogorov spectrum (Kolmogorov 1941) 
describing the spatial properties of the small-scale turbulence. The temporal behaviour of 
the small-scale turbulent dynamics has a frequency spectrum that is approximated by a 
Gaussian centred at zero frequency and with a width corresponding to the inverse of the 
correlation timescale of an eddy whose spatial extent (or wavenumber) is characterized by 
the mixing length. The correlation timescale is not a well-defined quantity, and therefore we 
scale a (well-defined) characteristic eddy turnover time with a correlation parameter A (cf. 
Balmforth 1992b), which we set to A = 1. 
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2.2. Mode peEik pcirameters 

In this section we describe how the resuhs of the stellar model computations were used 
to make predictions of the parameters of the p-mode peaks observed in the frequency power 
spectrum. For given values of a? and 6^ the model computations provided predictions of the 
linear damping rates, r], and the acoustic energy supply rates, P, of the radial p modes. The 
observed parameters of the mode peaks are formed from these quantities. The peak fwhm 
linewidths are given by 

A = 77/7r. (1) 
The mode velocity powers, V'^, are calculated from (e.g., Houdek et al. 1999): 

where / is the mode inertia. The maximum power spectral density, or height H, of a mode 
peak in the frequency power spectrum depends on the effective length of the dataset, T, and 
the hnewidth. A, via the formula (Fletcher et al. 2006): 

\ril ) TjT + 2 ttTA + 2' ^ ' 

When ttTA << 2 (i.e., when T << 2/7]) the mode is not resolved and power is confined 
largely in one bin of the frequency power spectrum, so that H . When, on the other 
hand, ttTA >> 2 (i.e., T >> 2/?]), the mode is well resolved and the Lorentzian shape of 
the peak may be inferred from the power in the bins occupied by the peak. For cases in 
the intermediate regime - where ttTA is neither much greater, or much smaller, than two - 
there is a gradual transition between H for the unresolved and fully resolved regimes. 

In other recent papers (e.g., Chaplin et al. 2005) we have gone on to use H to make com- 
parisons of theory with observation. In this paper we instead use data on heavily smoothed 
frequency power spectra, which give direct inference on (and therefore the amplitudes, 
V) as opposed to the H. However, we do still make use of the heights to calibrate the model 
computations, as we now go on to discuss. 



2.3. Notes on calibration 

We make three important points concerning the calibration. Firstly, the computations 
have been calibrated so that, for a model of the Sun, the average maximum power spectral 
density, H, of the five most prominent modes is the same as that observed in the real BiSON 
Sun-as-a-star data. By using an average over several modes, as opposed to taking the H of 
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just the strongest mode, we seek to stabilize the cahbration against small-scale fluctuations 
in the computations (Chaplin et al. 2007). Furthermore, by using data on H, as opposed to 
V"^ (or V), we maintain consistency with our previous work (e.g., see Chaplin et al. 2008). 
To summarize our first point: data on the Sun serve as a reference calibration for the model 
computations. 

Secondly, it is important to recognise that the reference calibration had to be performed 
independently for each of the non-local parameter choices a? = h'^ = 600 and o? = }p' = 300, 
respectively. This is because the changes to the values of the non-local parameters can have 
a significant impact on the computed damping rates, which in turn affects the absolute 
magnitudes of the predicted velocity powers, V"^ (Equation [2]) and therefore the heights, H 
(Equation [3]). 

Thirdly, when it comes to comparing the results of the pulsation computations with the 
observational data on a Cen B, we must remember that there will be instrument-dependent 
differences in the observed velocity amplitudes, due to the use of different spectral lines in 
the Doppler velocity observations. We have folded into the calibration of the model data 
the fact that the amphtudes of the solar p modes measured by the stellar techniques are a 
factor 1.07-times smaller than the amplitudes measured by BiSON (Kjeldsen et al., 2008). 



3. Results 

3.1. Data from the observations 

The left-hand panel of Fig. [1] shows the observed frequency power spectrum of a Cen B, 
plotted on a logarithmic scale. The spectrum was computed from a few nights of Doppler 
velocity data collected by Kjeldsen et al. (2005). The dark solid line is a smoothed spectrum 
given by applying to the raw spectrum a Gaussian filter of width 4At/, where is the 
large frequency spacing between consecutive overtones (here 162 /iHz). The dashed line is 
a smooth estimate of the background power spectral density. It was obtained by fitting, 
in regions outside the range occupied by the p modes, a second-order polynomial to the 
logarithm of power versus the logarithm of frequency. 

We used the Gaussian smoothed spectrum and the background fit to estimate the mode 
amplitudes, V . This was done by following the recipe outlined in Kjeldsen et al. (2008)0. 



^Note that we get mode amplitudes for a Cen B and the Sun that are a few percent lower than in 
Kjeldsen et al. These differences arise from differences in the fitting-function that was used to estimate the 
background. The estimation of the background is the largest source of uncertainty for the method. 
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In summary, we began by subtracting the background fit from tlie Gaussian-smootlied spec- 
trum. Tfie residuals tliereby obtained were converted to units of power per Hertz, multiplied 
by the large frequency spacing Au, and finally divided by a constant factor (see factors in 
Table 1 of Kjeldsen et al. 2008) to allow for the effective number of modes in each slice 
Au of the spectrum. This gave observational estimates of the radial mode amplitudes, Vobs, 
which are plotted as a thick grey line in each panel of Fig. [2l The thin grey lines are an 
estimate of the uncertainty envelope on the observed amplitudes. These uncertainties were 
estimated by using as a guide the results of analyzing many independent, short segments of 
BiSON Sun-as-a-star data, as will be explained below. 

We have used results on a parallel analysis of the BiSON Sun-as-a-star data as a belt- 
and-braces check on the procedures and results. The complete BiSON timeseries that we 
used is 4752 days long. This was split into independent segments of length 5 days (a length 
similar to the a Cen B timeseries), and the analysis procedures that were applied to the 
a Cen B frequency power spectrum were applied to the frequency power spectrum of each 
segment. The right-hand panel of Fig. [1] shows the frequency power spectrum of one of the 
5-day BiSON segments (the layout and linestyles are the same as in the left-hand panel, and 
the spectrum was again smoothed over 4Az/, but with Au = 135 /iHz for the Sun). 

The thick grey lines in both panels of Fig. [3] show the mean of the amplitudes that were 
obtained from the ~ 950 independent 5-day BiSON segments, using the method of Kjeldsen 
et al. outlined above. The thin grey lines mark plus and minus the rms of the amplitudes. 
Since the total epoch covered by the data spans more than one 11-yr cycle of solar activity, 
these uncertainties reflect variation of the estimated amplitudes from short-term stochastic 
variability and long-term solar cycle variability. There will also be a contribution from the 
finite signal-to-noise ratio of the observations. 

The fractional uncertainties in the estimated 5-day BiSON amplitudes are about 15% 
in the middle of the spectrum; and about 25 % at the extreme frequencies ~ 2000 and 
~ 5000 //Hz, respectively. In the light of these values, we have assumed the estimated a Cen 
B amplitudes are all determined to a fractional precision of 20 %. This value was used to 
make the uncertainty envelope for the estimated a Cen B amplitudes plotted in both panels 

of Fig. m 

With estimates of the uncertainties on Vobs now in hand, next we ask the question: how 
robust are the Vobs likely to be for a Cen B? We can obtain some insight by testing the 
robustness of the 5-day BiSON Sun-as-a-star amplitudes. This test may be accomplished 
by comparing the 5-day estimates with estimates from a "peak-bagging" analysis of the 
entire 4752-day BiSON timeseries. The BiSON mode amplitudes are usually estimated 
using the peak-bagging fitting techniques (e.g., see Chaplin et al. 2006). Peak-bagging 
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involves maximum-likelihood fitting of mode peaks in the frequency power spectrum to 
multi-parameter fitting models, where individual mode peaks are represented by Lorentzian- 
like functions. The points with error bars in Fig. [3] show estimated amplitudes Vbag from a 
full peak-bagging analysis of the 4752-day BiSON timeseries, in which the mode peaks were 
fitted in the high-resolution frequency power spectrum of the complete timeseries. Since these 
peak-bagging estimates are extremely precise, and are also expected to be fairly accurate, 
they serve to provide a robust cross-check of the 5-day estimates. As we can see, the 5- 
day BiSON estimates are in good agreement with the peak-bagging BiSON estimates. This 
comparison suggests we may also expect to have reasonable confidence in the observed a Cen 
B amplitudes. 



3.2. Comparison of theory with observation 

Predictions from the pulsation computations of the p-mode velocity amplitudes, V ^ are 
shown as dotted lines in Fig. [2] (a Cen B) and Fig. [3] (Sun). The predicted amplitudes plotted 
in the left-hand panels of the figures are for c? = }? = 600, while those in the right-hand 
panels are for = 1? = 300. 

First, let us compare the modelled and observed amplitudes for a Cen B (Fig. [2]). When 
a} = }? = 600 (left-hand panel) the model amplitudes form a very broad plateau at frequen- 
cies > 3700 /iHz, and the level of this plateau rises slowly with increasing frequency. While 
the match between the predicted amplitudes and the observed amplitudes is reasonable at 
frequencies < 4500 /zHz, this is demonstrably not so at higher frequencies, where the pre- 
dicted amplitudes are significantly higher than the observed amplitudes. When d} = }? = 300 
the predicted amplitudes give a better match, on average, to the observed amplitudes, but 
some overestimation at higher frequencies remains. 

High-frequency disagreement between theory and observation is also seen in Fig. [3] for 
the Sun. When a} = }? = 600, we again see a pronounced disagreement at high frequencies 
between the predicted amplitudes and the observed amplitudes. When = = 300 the 
level of disagreement is less pronounced, but is nevertheless still present. 

Let us summarize the main points from Figs. [2] and [H when o? = = 600, theoret- 
ically computed high-frequency velocity amplitudes for a Cen B and the Sun significantly 
overestimate the observed velocity amplitudes. Although not as severe, this overestimation 
persists at = 6^ = 300. What might be the cause of this disagreement between theory 
and observation? That is the question we turn to next. 
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4. Inference on the p-mode linewidths 

Provided we trust the mode amphtudes Vobs estimated from the a Cen B frequency 
power spectrum - and the 5-day BiSON analysis above suggests the estimates should be 
reasonable - there are two possible ways out of the problem. Equation [2] implies that: 

5V 6P 67] 

— ~ -. (4) 

V 2P 2r] ^ ^ 

So, either the model damping rates t] - and therefore the model linewidths A - are too weak 
at high frequencies to explain the observed Vobs, or the modelled acoustic supply rates, P, 
are too strong (or there is a combination of the two effects). 

There is good evidence from analysis of Sun-as-a-star data (e.g., Chaplin et al. 2005; 
Houdek 2006) that a significant part of the disagreement may come from problems computing 
T) (and therefore A). We can double-check the validity of this assumption here for the Sun, 
because we have precise BiSON data on the linewidths and powers of the solar p modes from 
the peak-bagging analysis. These peak-bagging data may be used as a precise and accurate 
reference against which to check the quality of the results obtained on the short 5-day BiSON 
segments. 



4.1. The Sun-as-a-star linev^idths 

We base our solar check on Equation [2J Again, it tells us that cx P/A. If we take 
the ratio of the observed to the model-computed velocity powers, we will therefore have: 



obs \ / -Tobs 



VP 



^obs 



(5) 



In the above, all observed quantities for the star, which we assume come from analysis of 
the short 5-day segments, carry the suffix "obs"; and the model-predicted quantities are 
suffix-free. Rearrangement of the above gives: 



Aob. = ( ^ ) ( ) A. (6) 



Let us suppose that all the problems in the modelling lie in computation of the linewidths, 
A. The implication is then, of course, that the model- predicted values of the acoustic supply 
rates, P, are accurate, i.e., P = Fobs- Equation would then simplify to 



^obs 



^obs 



A. (7) 
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How well do the solar model computations and 5-day observations match Equation [TJ i.e., 
how well is the relation P = Pohs satisfied? Fig. H] shows plots of the inferred linewidths 
Aobs of the Sun (thick grey line). The Aobs were computed from the observed 5-day Vobs 
in Section 13.11 above and the model predictions of V"^ and A. The left-hand panel shows 
results for when the model computations use a? = h"^ = 600; while the right-hand panel 
shows results with a? = 1)^ = 300. The thin grey lines show the uncertainty envelopes on the 
Aobs, which come from the uncertainties on Vobs shown in Fig. [3l 

Also plotted as points with associated best-fitting uncertainties in Fig. H] are the fitted 
linewidths Abag from a full 4752-day peak-bagging analysis of the BiSON Sun-as-a-star data. 
For our purposes here, we may regard the peak-bagging linewidths as being good measures 
of the true linewidths of the Sun. Our check on whether P = Fobs therefore amounts to 
seeing if the 5-day Aobs are a good match to the Abag- This does assume that there are no 
significant biases in the estimated Vobs, as was suggested by the good agreement of the Vobs 
and Vbag in Fig. [31 

At a first glance, the most striking aspect of the solar plots is indeed the encouraging 
level of agreement between the Aobs and the Abag- Fig. O shows the comparison in more 
detail. Here, we have plotted the fractional differences between the 5-day linewidths and 
the peak-bagging linewidths. When a? = }p' = 600, the Aobs over the main part of the solar 
p-mode spectrum are seen to be about 10 % higher on average than the Abag, while they are 
about 35 % lower when = b'^ = 300. The implication is that in both cases there is actually 
an offset between P and Fobs- However, what we can say is that changes in the differences as 
a function of frequency are not significant over the main part of the p-mode spectrum, given 
the observational uncertainties (we can disregard the differences at the lowest frequencies, 
where the error bars are very large, and the 5-day spectra have insufficient resolution to 
give robust estimates of Aobs in this part of the spectrum). We are therefore in a position 
to conclude the following: the shape in frequency of the Aobs is a reasonable match to the 
shape of the Abag at the level of precision of 5-day timeseries. 

The above suggests we may use Equation [7] to infer the variation of the linewidths as 
a function of frequency, without the need to measure those linewidths directly. There will 
remain some uncertainty over the absolute calibration of the linewidths. Let us now apply 
Equation [7] to the data on aCen B. 
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4.2. Inference on linewidths of a Cen B 

We have used the observed estimates of Vobs from the a Cen B spectrum, together 
with the theoretical computations of V and A for the star, to give the inferred hnewidths 
Aobs plotted in Fig. [61 We do not require an explicit estimate of Fobs, because we rely on 
the assumption (verified above for the Sun) that the theoretical P has the same shape in 
frequency as Fobs- 

There is no obvious reason why we should expect the acoustic supply-rate computations, 
P, to be any less valid for a Cen B than they are for the Sun (this might not have been 
so had a Cen B been somewhat hotter than the Sun, or much cooler than it actually is; 
see Houdek 2006). At the very least we may therefore have reasonable confidence in the 
shapes of the inferred linewidth curves for the star, assuming any deviations from a constant 
scaling offset are no more severe than they are for the Sun. Both curves in Fig. [6] appear 
to show a plateau in the linewidths at 4000 /zHz, like that seen in the solar linewidth curve. 
Other obvious features of the curves - a decrease of the linewidths at lower frequencies, 
and an increase at higher frequencies - are also Sun-like in nature. It is worth stressing 
that the linewidth curves have very similar shapes at = \? = 600 (left-hand panel) and 
= }? = 300 (right-hand panel). 

We have also shown on Fig. E] observational linewidth estimates, which were obtained 
by Kjeldsen et al. (2005) in two parts of the p-mode spectrum (points with associated uncer- 
tainties). These estimates were obtained from the same data but in a different way, namely 
by measuring the scatter of p-mode frequencies around smooth ridges in the echelle diagram. 
Comparison of the inferred Aobs with the Kjeldsen et al. estimates allows us to place con- 
straints on the absolute values of the linewidths for a Cen B. Such a comparison suggests 
there is actually little to choose between the Aobs and the Kjeldsen et al. estimates when the 
energy supply rates are computed using a} = \? = 600 (giving the inferred linewidths in the 
left-hand panel of Fig. [6]) or = 6^ = 300 (giving the inferred linewidths in the right-hand 
panel). However, the a?' = \? = 300 result does give a slightly better match. Inspection 
of the Aobs linewidth curve then implies that the a Cen B linewidths take values of about 
0.4 /iHz at a frequency of 3500 /xHz, about 0.8 /xHz at a frequency of 4000 /xHz, and about 
2.5 yuHz at a frequency of 4500 yuHz. These linewidths turn out to agree reasonably well with 
the solar linewidths, if the mode frequencies are multiplied by the ratio of the acoustic cut-off 
frequencies of the two stars (which in effect scales the a Cen B frequencies down to those 
shown by the Sun). In short, the a Cen B linewidth curve is similar to a frequency-scaled 
version of its solar cousin. 
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5. Discussion 

In this paper we presented an analysis of tlie amplitudes and linewidths of the low-degree, 
Sun-like p modes displayed by the star a Cen B. These data were extracted from only a few 
nights of Doppler velocity observations collected on the star by Kjeldsen et al. (2005), and 
were compared with theoretical predictions of the stochastic excitation and damping rates of 
the p modes. We also performed a parallel analysis of Sun-as-a-star p-mode data collected 
by the ground-based BiSON. The very long BiSON timeseries allowed us to validate the 
analysis techniques. 

For the Sun, we found that model predictions of the mode amplitudes were significantly 
larger than the observed amplitudes in the high-frequency part of the p-mode spectrum. We 
were able to confirm that most of the disagreements for the Sun, which set in at frequencies 
> 3300 /iHz, are due to problems computing the damping rates, not the excitation rates, 
of the modes. The computed damping rates must be increased to explain the observed 
amplitudes. 

Similar disagreements are seen for a Cen B. Here, we do not have the luxury of cross- 
checking the analysis with results from very long datasets, since the latter at present do not 
exist [although we may hope to obtain such datasets in the future from the likes of SONG 
(Grundahl et al. 2007) and SIAMOIS (Mosser et al. 2007)]. However, by assuming that 
the disagreements for a Cen B are also due largely to problems with the model-predicted 
damping rates, as was shown to be the case for the Sun, we were able to demonstrate that 
the model linewidths must also be increased significantly in order to explain the observed 
amplitudes. The problems for a Cen B set in at frequencies > 4500 /^Hz. 

The conclusions above bear on a prediction of the pulsation computations mentioned 
in the Introduction (Section [1]): that stars cooler than about 5400 K will show a broad 
plateau or double-hump of power in their p-mode spectra. We showed that a Cen B does 
not have the predicted broad plateau in its observed p-mode spectrum. If the high-frequency 
damping rates are increased - as they must be to resolve the disagreement between theory 
and observation - the broad plateau or second high-frequency hump disappears in the cooler 
models. 

Finally, we also showed how, by making suitable comparisons of theory and observation, 
it is possible to make inference on the variation with frequency of the p-mode linewidths 
without directly measuring those linewidths, even if the observations come from only a few 
days of data. Use of independent measures of the a Cen B linewidths, made in two parts 
of its spectrum by Kjeldsen et al. (2005), allowed us to calibrate our inferred linewidth 
curve for a Cen B. We found that the resulting, calibrated linewidth curve is similar to a 



-14- 



frequency-scaled version of its solar cousin, with the scahng factor equal to the ratio of the 
respective acoustic cut-off frequencies of the two stars. The ratio of the frequencies at which 
the onset of high-frequency problems is seen in both stars is also given approximately by the 
same scaling factor. 
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Fig. 1. — Left-hand panel: Observed frequency power spectrum of ctCen B (Kjeldsen et al. 
2005). The dark sohd hne is the smoothed spectrum, while the dashed line is an estimate of 
the background. Right-hand panel: Observed frequency power spectrum of a 5-day segment 
of BiSON Sun-as-a-star data (linestyles as per left-hand panel). 
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Fig. 2. — Observed mode velocity amplitudes Vobs (thick grey line) and theoretically com- 
puted amphtudes V (dotted line) of a Cen B. The thin grey hues denote the estimated 
uncertainty envelope on the amplitudes (see text). Left-hand panel: theoretical predictions 
for — — 600. Right-hand panel: theoretical predictions for — — 300. 
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Fig. 3. — Observed mode velocity amplitudes Vobs (thick grey line) and theoretically com- 
puted amplitudes V (dotted line) of the Sun. The plotted amplitudes are the mean ampli- 
tudes from analyzing 950 independent 5-day segments of BiSON Sun-as-a-star data. The 
thin grey lines denote the uncertainty envelope on the amplitudes, from the scatter in the 
results from the 5-day segments. The symbols with error bars are amplitudes Vbag given by a 
"peak-bagging" analysis of the full 4752-day BiSON timeseries. Left-hand panel: theoretical 
predictions for = 6^ = 600. Right-hand panel: theoretical predictions for a? — h"^ — 300. 
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Fig. 4. — Inferred linewidths Aobs (thick grey line) and theoretically computed linewidths A 
(dotted line) of the Sun. The theoretically computed linewidths have been smoothed with 
a median smoothing filter of three consecutive values, as in Houdek & Gough (2002). The 
plotted inferred linewidths are the mean inferred linewidths from analyzing 950 independent 
5-day segments of BiSON Sun-as-a-star data. The thin grey hues denote the uncertainty 
envelope on the inferred linewidths, from the scatter in the results from the 5-day segments. 
The symbols with error bars are linewidths Abag given by a "peak-bagging" analysis of the 
full 4752-day BiSON timeseries. Left-hand panel: theoretical predictions for — — 600. 
Right-hand panel: theoretical predictions for = = 300. 




Fig. 5. — Results on the solar p-mode linewidths. Plotted are differences in the natural 
logarithms (i.e., absolute fractional differences) of the inferred linewidths, Aobs (which use 
data from the 5-day BiSON segments) and pealc-bagging linewidths, A^ag (from the frequency 
power spectrum of the full 4752-day BiSON timeseries). 
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Fig. 6. — Inferred linewidths Aobs (thick grey line) and theoretically computed linewidths A 
(dotted line) of a. Cen B. The theoretically computed linewidths have been smoothed with 
a median smoothing filter of three consecutive values, as in Houdek & Gough (2002). The 
thin grey lines denote the estimated uncertainty envelope on the inferred linewidths (see 
text). The symbols with error bars show observational linewidth estimates from Kjeldsen et 
al. (2005). Left-hand panel: theoretical predictions for — }? — 600. Right-hand panel: 
theoretical predictions for = \? — 300. 



